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(54) Producing medical images corrected for patient movement 

(57) An x-ray CT scanner acquires projection data from a series of views during a scan of a patient's chest Movement of 
the patient's chest due to respiration is also sensed (e.g. using an ultrasonic range finder 44, Fig 5) during the scan and this 
acquired motion data is employed along with a geometric model of chest motion to calculate factors which correct the 
acquired projection data and reduce motion artifacts in an image produced by back projecting the acquired projection data 
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PRODUCING MEDICAL IMAGES 

The present invention relates to medical imaging, and 
particularly, to the reduction of motion artifacts in images 
produced using a projection method of reconstruction. 
S There are a number of modalites used to produce medicai 

images. These include x-ray computed tomography (CT) , 
magnetic resonance imaging (MRI) , single photon emission 
computed tomography (SPECT) , and positron emission computed 
tomography (PET) methods. In all cases, the data used to 

10 reconstruct the desired image is acquired over a period of 

time in a scan comprised of a series of views. Each view is 
a snapshot of the patient from a different angle, or 
perspective, and a scan typically includes tens, or hundreds 
of views. In the case of x-ray CT the entire data set may be 

15 acquired in a few seconds, whereas an MRI scan typically 

requires a few minutes to complete. The methods used to 
reconstruct an image from such data sets presume that the 
patient is motionless during the entire scan and that the 
same fixed object is the subject of all acquired views. To 

20 the extent this is not true, artifacts such as ghosts, 

smearing and fuzziness appear in the reconstructed image. 

Efforts to reduce patient motion during a scan can 
significantly improve image quality. However, artifacts 
caused by respiration is a significant problem in chest scans 

25 where suspension of breathing is not possible or poor 

instructions are provided to the patient by the scanner 
operator. Children and comatose patients are routinely 
scanned with no attempt to synchronize respiration with 
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scanning, and it is expected in such cases that a number of 
poor quality images will be produced and will be discarded. 

One approach to reducing motion artifacts in medical 
images is to retrospectively correct the acquired data to 
offset the affects of motion. One such method, for example, 
is disclosed in U.S. Patent No. 4, 937,526 and is applied to 
acquired MRI data. The corrections that are made may be 
determined from an examination of the acquired raw data 
itself, or additional information, such as a signal from a 
cardiac monitor or a respiration monitor, may be used. The 
manner in which the corrections are made to the acquired raw 
data is determined by the particular reconstruction technique 
that is used. In the above patent, for example, a 2D Fourier 
transformation is used to reconstruct an image from the 
acquired MRI data, and the correction methods disclosed are 
limited to that technique. 

The back projection method for image reconstruction is 
employed to some extent in all computed medical imaging 
modalites. It is the predominant method used in x-ray CT, 
and there is a need to correct acquired data used in 
projection imaging for the effects of patient motion. 

Summary of t_h* Tnvpnti^ 

The present invention relates to a method for correcting 
data used to reconstruct medical images using a back 
projection technique, and particularly, to a method for 
reducing artifacts in the reconstructed image due to 
respiratory motion of the patient during the data 
acquisition. More specifically, the method includes 
acquiring the image data during a scan comprised of a series 
of views, acquiring motion data with each view that measures 
a parameter associated with patient respiration, correcting 
the acquired image data as a function of the acquired motion 
data and a geometric model of the patient's chest cavity 



-3- 

during respiration, and reconstructing an image by back 
projecting the corrected acquired image data, 

A general object of the invention is to reduce artifacts 
produced in back projection images caused by patient 
5 respiration. A model of how the patient's chest changes in 

size and shape as a function of a measured parameter has been 
developed. By acquiring the measured parameter along with 
each view in the scan, the size and shape of the patient's 
chest as seen by each view can be determined. During the 

10 reconstruction/ the acquired data is first corrected by a 

weighting factor and is then used in a back projection 
process which "sees" the patient's chest in a single 
reference position. 

A more specific object of the invention is to provide a 

15 method for weighting image data acquired by an x-ray CT and 
reconstructing an image from the weighted image data using a 
back projection technique. The back projection process is 
altered by a displacement factor to properly reconstruct the 
image and the specific weighting factors and displacement 

20 factors for both parallel beam and fan-beam CT scanners have 

been developed. 

The foregoing and other objects and advantages of the 
invention will appear from the following description. In the 
description, reference is made to the accompanying drawings 

25 which form a part hereof, and in which there is shown by way 

of illustration a preferred embodiment of the invention. 
Such embodiment does not necessarily represent the full scope 
of the invention, however, and reference is made therefore to 
the claims herein for interpreting the scope of the 

30 invention. 

Brief nftscription of the Drawings 

Fig. 1 is a perspective view of an x-ray CT scanner and 

patient; 
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Figs. 2A and 2B are schematic drawings of a parallel beam and 
a fan-beam scanning assembly on the scanner of Fig. i ; 

Fig. 3 is an electrical block diagram of the scanner of Fig. 
l; 

5 Fig. 4 is a graphic representation of one view, of data 
acquired by the scanner of Fig. 1; 
Fig. 5 is a schematic representation of a patient's chest 
cavity showing the orientation of a range finder that 
forms part of the scanner of Fig. 1; 
10 Fig. 6 is a flow chart of the program executed by the CT 

scanner of Fig. 1 to carry out the preferred embodiment 
; of the invention; and 

Fig. 7 is a schematic representation of one of the steps 
employed by the program of Fig. 6. 
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Description of the p^^t-^h E mhnHimon*- 

While the present invention may be applied to many 
different imaging systems that employ back projection image 
reconstruction methods, the preferred embodiment is employed 
in an x-ray CT scanner such as that illustrated in Fig. 1. 

As shown in Fig. 1, a CT scanner used to produce images 
of the human anatomy includes a patient table 10 which can be 
positioned within the aperture 11 of a gantry 12. A source 
of highly columinated x-rays 13 is mounted within the gantry 
12 to one side of its aperture 11, and one or more detectors 
14 are mounted to the other side of the aperture. The x-ray 
source 13 and detectors 14 are revolved about the aperture 11 
during a scan of the patient to obtain x-ray attenuation 
measurements from many different angles. 

A complete scan of the patient is comprised of a set of 
x-ray attenuation measurements which are made at discrete 
angular orientations of the x-ray source 13 and detector 14. 
Each such set of measurements is referred to in the art as a 
"view" and the results of each such set of measurements is a 
transmission profile, or projection. As shown in Fig. 2A, 
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the set of measurements in each view may be obtained by 
simultaneously translating the x-ray source 13 and detector 
14 across the acquisition field of view, as indicated by 
arrows 15. As the devices 13 and 14 are translated, a series 
5 of x-ray attenuation measurements are made through the 

patient and the resulting set of data provides a transmission 
profile at one angular orientation (9) . The angular 

orientation of the devises 13 and 14 is then changed (for 
example, 1') and another view is acquired. These are known 
10 in the art as parallel beam projections. An alternative 

structure for acquiring each transmission profile is shown in 
Fig. 2B. In this construction, the x-ray source 13 produces 
a fan-shaped beam which passes through the patient and 
impinges on an array of detectors 14, The detectors 14 can 
15 be curved as shown in Fig. 2B, or they can be aligned in a 

straight line (not shown in the drawings). Each detector 14 
in this array produces a separate attenuation signal and the 
signals from all the detectors 14 are separately acquired to 
produce the transmission profile for the indicated angular 
20 orientation. As in the first structure, the x-ray source 13 

and detector array 14 are then rotated to a different angular 
orientation and the next transmission profile is acquired. 

As the data is acquired for each transmission profile, 
the signals are filtered, corrected and digitized for storage 
25 in a computer memory. These steps are referred to in the art 

collectively as preprocessing" and they can be performed in 
real time as the data is being acquired. The acquired 
transmission profiles are then used to reconstruct an image 
which indicates the x-ray attenuation coefficient of each 
30 voxel in the reconstruction field of view. These attenuation 
coefficients are converted to integers called r, CT numbers", 
which are used to control the brightness of a corresponding 
pixel on a CRT display. An image which reveals the 
anatomical structures in a slice taken through the patient is 
35 thus produced. 
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Referring particularly to Fig. 3, the operation of the 
CT system is controlled by a programmed data processing 
system 25 which includes a computer processor 26 and a disc 
memory 27. The disc memory 27 stores the programs the 
computer processor 26 uses in patient scanning and in image 
reconstruction and display. It also stores on a short-term 
basis the acquired data and the reconstructed image data. 
The computer processor includes a general purpose 
minicomputer with input and output ports suitable for 
connection to the other system elements as shown. It also 
includes an array processor such as that disclosed in U.S. 
Patent No. 4,4 94,141 which is encorporated herein by 
reference . 

An output port on the computer processor 26 connects to 
an x-ray control circuit 28, which in turn controls the x-ray 
tube 13. The high voltage on the x-ray tube 13 is controlled 
and its cathode current is controlled to provide the correct 
dosage. The high voltage and cathode current are selected by 
an operator who enters the desired values through an operator 
console 30 and the computer processor 26 directs the 
production of the x-rays in accordance with its scan program. 

The x-rays are dispersed in a fan-shape as described 
above and received by the array of detectors 14 mounted on 
the opposite side of the gantry aperture 11. Each individual 
25 cell, or detector element, examines a single ray originating 
from the x-ray tube 13 and traversing a straight line path 
through a patient located in the aperture 11. The currents 
formed in each detector element are collected as an analog 
electrical signal and converted into a digital number by A/D 
30 converters in a data acquisition system 31. The digitized 

measurements from all the detectors is a complete view. U.S. 
Patent Nos. 4,112,303 and 4,115,695 disclose details of the 
gantry construction, U.S. Patent No. 4,707,607 discloses the 
details of the detector array 14, and the data acquisition 
35 system is disclosed in U.S. Patent No. 4,583,240. All of 
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these patents are incorporated herein by reference. The 
digitized signals are input to the computer processor 26. 

The digitized attenuation measurements from the data 
acquisition system 31 are preprocessed in a well-known manner 
to compensate for "dark currents", for uneven detector cell 
sensitivities and gains, and for variations in x-ray beam 
intensity throughout the scan. This is followed by beam 
hardening corrections and conversion of the data to 
logarithmic form so that each measured value represents a 
line integral of the x-ray beam attenuation. This 
preprocessing is performed in real time as the scan is being 
conducted, and as shown in Fig. 4 f each view is comprised of 
a set of attenuation values 32 which define a transmission 
profile, or projection, indicated by dashed line 34. 

In addition to the transmission profile data 34, two 
other pieces of information are input during the acquisition 
of each view. The first is the angle (G) which indicates the 

angular orientation of the x-ray source 13 and detectors 14 
with respect to the vertical reference axis* Typically, for 
example, the views are acquired at 1* increments over a range 
of 180*. 

The second piece of information acquired with each view 
is a distance value (D) that is indicative of the position of 
the patient's chest cavity and which is a parameter employed 
in a geometric model of the chest cavity during respiration. 
Referring particularly to Fig. 5, this geometric model is 
illustrated by a schematic cross-section taken in a 
transverse plane through the patient's chest as indicated at 
40. As the patient breathes the posterior abdominal wall 41 
which rests on the supporting table 10 does not move any 
significant amount, whereas the anterior abdominal wall 42 
moves vertically as indicated by the dashed line 43. As will 
be described below, the size and shape of the patient's chest 
cavity at any point in the respiratory cycle can be 
approximated by monitoring the vertical position of the 
anterior wall 42. Accordingly, an ultrasonic range finder 44 
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is mounted to the gantry 12 and is positioned to measure the 
vertical distance (D) to the patient's ohest . This measured 
(D) is input to the computer processor 26 (Fig. 3) along with 
each view of acquired data. 
5 Referring still to Fig. 5, as the patient breathes and 

the anterior abdominal wall 42 moves up and down, the 
contents of the chest cavity magnify and shrink along the 
vertical axis (y) . This magnification does not occur about 
the center near y = 0, but instead, at the posterior wall 41 
10 located at y - -y p . There is virtually no magnification 

along the horizontal (x) , but to the extent that there is, it 
occurs about the center of the chest cavity at x=0 (if the 
patient is centered on the table 10) . Using this model and 
the measured parameter D, the present invention corrects the 
acquired projection data and employs the corrected projection 
data in the back projection reconstruction process such that 
the chest cavity appears stationary in a reference position 
during the entire scan. As a result, motion artifacts are 
significantly reduced or eliminated. 

The corrections to the acquired data and the manner in 
which the corrected data is employed in the back projection 
process has been determined for a general case in which 
patient motion produces magnification along two axes (x and 
y) and the point about which magnification occurs is shifted, 
or offset, from the origin (x - 0, y - 0) . The correction 
factors for this generalized case have been determined and 
will now be described. 

Let f (x,y) be the cross-section of the patient which is 
to be reconstructed. A magnified and shifted version of this 
same cross-section during various stages of the respiratory 
cycle is given by: 
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f(x,y) - flax + p x x, <Xy + (Jyy) 

where p x and py are magnification factors along the respect: 
x and y axes and a x and Oy are shift factors. For parallel 



[1] 
ive 



beam projection data the formula for projection of this 
magnified image at gantry rotational angle 9 is given by: 



p'(6,t) = j j f'(x,y) 8(t-x cos0-y sin6)-dxdy [2 



— OO —CO 



where 6(t) is the Dirac-delta function known to those skille 

in this art. The Fourier transform of this projection can L 
found: 

S'(9,co) = °j f <a x +P x x,a y +p y y)e -12m (x cosG+y sin6) dx dy [3 

— oo 

Now make the following changes of variables 
x- - Ct x+ p x x 

y' - a y +p y y S 
When (4] is used in [3], the following is obtained 

„ - j2«co| — cosG+^-sine j2Jld ~ cosG+^inO 

S'(e,(0) = f <XfY) e >-Px Py J e LPx Py J dx dy 

-i My 

Let F(u,v) be the two-dimensional Fourier transform of 
f(x,y). Then it is seen that 

j2ntd — cose+^3in9 
e LPx Py J cos0 ^ sin el 

s ' ,e ' a ' pa — F L"ftT'^rJ 16 
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z\:t:tj: a version ° f the ^ — ia 

^z^rz^:^z:^ ifi - and — 

p r o jecti0 n at gantry po -r e t r:;::i :rr ™ of the 

dimensional Fourier transform of the object function at angle 



9- - tan " 1 H* tan6 



P: 



after a phase term and a scaling factor have been removed 
Equations [6] and [7] could form the basis of a 

IIZZT^ * ^^nsional 

If 2 c T m ° f WhiCh r6SUltS f " m -pping 

dat L e r Si0nal tS -*« ° f .projection 

the T ! tranSf ° rm ° f the Patient - is not 

the preferred method of reconstruction, mstead, a filtered 
bac* projection method has been developed and will now be 
described in detail. 

now be d^f/ 30 " Pr ° jeCti ° n construction, formula will 
be derived for reconstructing f (x , yJ using p , ^ t) _ The 

inverse Fourier transform of F(u,v>, i s given by 



f (x,y> - j j F <u,v)e3 2 *<*« + vy) dudv 



[8] 



Consider the following change of variables 



U " t C ° S9 [91 



CO 

v = — - sine 



-11- 



The components of the Jacobian in this change of variables 
are given by 



du 
3d) 

9u 
89 

9v 
dco 

dv 

ae 



cos8 
Px 

CO sin9 
Px 

sin8 
co cos6 



cop 



• cos9 

a 



Px 2 



Pv 



coP' v sin8 
Pv 2 



where p f x and P' y are the derivatives of p x and py wi 
respect to 0. The Jacobian, J(u,v,co,6) , can be dete: 
using [10] resulting in 



th 

rmined 
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J(u,v,co,e) - Icolg (8) 



r i i - 



where 



g(9) - 



l + 



sin29 



Px 



[12; 



The value g(9) is a weighting factor which is applied to 
the projection data at gantry position 9. In the parallel 
15 beam acquisition/ this weighting factor is a constant value 

which is applied to each attenuation value 32 in the 
projection profile 34 at the position 0 (Fig. 4). In some 

situations where the density of the object being imaged does 
not change as it is magnified, this weighting factor is 
20 modified to more closely approximate the geometric model of 

motion in equation [1] . More specifically, in such cases the 
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weighting factor g(6) should be multiplied by the value p y p x 
The derivatives P' x and P' y can be calculated with numerical 
differences using the adjacent values of p x and 0y. 



Using [9 J and [11], [8] becomes 



5 fix vl=," T r I" " cos8 " 3ine1 , , a j H f "^f" 3 " 9 ] 

5 f(x,y) =j J F [ « jj~ J la>|g(8) e LPx J dfflde 



When [6] is used, [13] reduces to 



d9 [14] 



where 



q e (t) = JS'{6,0)) lco|g(0)ei 2,lCDt dO) [15: 



Equation [14] represents a filtered back projection formula 
for reconstruction of parallel projections that are acquired 
from the magnified and shifted object. The formula is valid 
for reconstructing any point in the x-y plane. 

The image f(x,y) can be reconstructed, therefore, from a 
15 set of parallel beam projections acquired over a range of 

180' gantry positions by modifying the conventional data 
acquisition and reconstruction method in the following 
manner. First, for each projection of the patient, not only 
is the attenuation profile data acquired, but also, the 
parameters <X X , p x , Oy and p y are measured. In the preferred 
embodiment, is set to zero and p x is set to one since the 
patient is usually centered on the table 10 and there is very 
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little magnification of the chest cavity along the x axis 
during respiration. Only (iy and <Xy are required, therefore, 

to significantly reduce motion artifacts in the chest cavity 
and these are measured indirectly. As shown in Fig. 5, the 
distance between y = 0 and the point about which 
magnification occurs is fixed at -y p . Also, the distance 
between the range finder 4 4 arid this same point is fixed. As 
a result, the values for py and Oy can be calculated from 

these fixed values and the measurement (D) produced by the 
range finder 44 as follows: 

p y - i/d + (D - D 0 )/(yo + y P )> us: 

OCy = -y p (1 - Py) [17! 

where y 0 is a reference position for the anterior chest wail 
which is selectable by the operator and which determines the 
shape and si2e of the reconstructed image, and Do is the 
range. finder measurement at this reference position* 
Consequently, the distance measurement (D) is acquired alcr.g 
with each projection profile and this measured parameter is 
sufficient to indicate the shape and size of the patient's 
chest cavity at the moment the view was acquired. 

The preferred embodiment of the invention will now be 
described with reference to the flow chart of Fig. 6. While 
"most of the steps are carried out in dedicated hardware so 
that the processing can be carried out in "real time", the 
process itself is controlled by a program executed by the 
computer 26 which performs the scan. 

This control program is entered at 75 and the CT system 
is initialized at process block 76 to acquire the data for 
the first view. This includes receiving input data from the 
operator such as the reference chest position y Q and 
reference range finder distance Do, and orientation of the 
gantry to the desired starting position of 8 - -90-. A loop 
is then entered in which the profile data for the first view 
is acquired and preprocessed as indicated at block 77. The 
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^stance measurement <D> from the range finder 44 (Fig. 5, - s 
acquired at process block 78 and the values for Py anc I J a ~ 
calculated at process block 79 using the equations (16] and " 
[17]. The acquired projection data is filtered in the usua'- 
fashion at process block 80 and then it is weighted at 
process block 81 by multiplying each value in the profit 
data set by the weighting factor g(9 ) calculated in " 
accordance with equation [12]. 

The corrected projection data is employed in 
reconstructing an image using the back projection technics 
as indicated at process block 82. As indicated by equation 
[14 J, however, this process is modified to account for 
motion. Referring particularly to Fig. 7, a 512 by 512 pr- 
image 85 is created by determining which of the values in f- 
corrected and filtered projection data set 86 contribute -0 
the brightness value of the pixel located at (x,y) . m a 
parallel beam acquisition the conventional back projection 
formula for determining which value (t) to use is as follows: 

t - x cosG + y sinQ 15 - 

where <x,y> is the location of the pixel, 6 is the projection 
angle for the view, and t is the location in the projection 
data set from which an attenuation value 87 is read. Th<s 
standard back projection is shown in Fig. 7 by the dashed 
line 88. Typically, t is located between two samples in t-e 
acquired data set and interpolation is used to determine a 
more accurate value to be added to the CT number for pixel 
<*,y> • For each projection, all of the pixels in the image 
85 are processed in this fashion to determine the 
contribution to their accumulated CT numbers. 

To practice the present invention this back projection 
technique is changed to select a different value <f> from 
the corrected projection data set 86. This selection is made 
as follows: 

t 1 * x' cos6 + y' sin8 r:3' 
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where 

x' - (x - a K ) / p x [2C3 

y - (y - ay) / p y [21: 

That: is, the back projection process is modified by a 
5 displacement factor that is determined by the values a x , a y , 
Px and Py. This change is illustrated in Fig. 7 where (x,y) 
is the pixel in the reference image being reconstructed, 
(x'^y'l is the location of the same point in the patient, at 
the time the projection data was actually acquired, and the 

10 attenuation value 89 is the value selected by equation [19] . 

In other words, the geometric model and the motion parameter 
D indicate that the attenuation value to be used at the pixel 
(x,y) from the projection data 86 is the value 89 at t' 
rather than the attenuation value 87. After the contribution 

15 to each pixel in the image has been computed, the system 

loops at decision block 90 to advance the gantry and acquire 
and process that data for the next view. When 180* of data 
have been acquired and processed in this manner, the scan is 
complete and the image data 85 is displayed at process block 

20 91. The CT numbers in the image data array 85 are scaled and 
processed in the normal fashion to produce an image of the 
desired brightness level and range. 

The teaching of the present invention is also applicable 
to fan-beam CT scanners which employ the back projection 

25 technique of image reconstruction. As in the preferred 

embodiment described above, the acquired projection data is 
corrected by a weighting factor that is determined by the 
magnification and shift values Ox, Oy, p x and Py, and the back 

projection process is modified by a displacement factor that 
30 is also determined by these same values. The calculation of 
the weighting factors and displacement factors are different, 
however, due to the different geometry. 
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For a flat detector fan-beam reconstruction the 
weighting factor used in process block 81 in Fig. 6 is as 
follows : 



G(S,d) = 



1 + 



sin 2 (a+tan" 1 S/R) 



R2 
2 PxPv 



[21} 



15 



5 where a is the rotational position of the gantry, s is the 

position of the x-ray detector which is being weighted with 
respect to the center of the detector array, R is the 
distance between the x-ray source and the center of the 
detector array, and 

10 z - V R2 + s 2 . 

In contrast to the parallel beam acquisition, this weighting 
factor not only varies as a function of gantry position a, 
but also as a function of the location of the detector in the 
flat array. During the back projection process of block 32 a 
different formula than equation [18] is used for selecting 
the proper attenuation value for each pixel (x,y) . Many back 
projection formulas are known in the art such as that 
described in U.S. Patent No. 4,812,983 entitled "Method and 
Means of Correcting For a Shift in the Center of Rotation of 
a Rotating Fan-Beam CT System- which is encorporated herein 
by reference. Regardless of the formula used, the 
displacement factor of the present invention is applied by 
substituting the values of x- and y given above in equations 
[20] and (21] for the values of x and y respectively in the 
25 particular back projection formula used. 

For a curved detector fan-beam reconstruction the 
weighting factor used in process block 81 in Fig. 6 is as 
follows : 



20 
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G(Ot,Y> = 



1 + 



sin 2 (OC+y) 



R cosy 



[22] 



where a is the rotational position of the gantry, R is the • 
distance between the x-ray source and the central axis of 
rotation of the gantry, and y is the angle as measured at the 
5 x-ray source between the central array detector and the 

detector whose signal is being weighted. During the back 
projection process of block 82, the values of x' and y' giver, 
above in equations [20] and [21] are substituted for the 
values of x and y respectively in the formula used for back 

10 projection. 

While the theory indicates that the weighting factors 
must be applied to correct the projection data before it is 
used to reconstruct an image according to the present 
invention, experimental results have shown that this is not 

15 always required. In many cases, a substantial reduction in 

motion artifacts can be achieved without applying the 
weighting factor and only applying the displacement factor to 
the back projection process. 

It should be apparent to those skilled in the art that 

20 the present invention is applicable to many different back 

projection reconstruction techniques. Regardless of the back 
projection technique used, a weighting factor can be 
calculated for each acquired attenuation value in the data 
set and the back projection process can be modified by 

25 substituting the displacement factors of equations [20] and 

[21] into the back projection formula. This is true 
regardless of the modality used to acquire the projection 
data. Thus, for example, projection data acquired with PET, 
MRI or SPECT scanners can be corrected for patient motion 

30 according to the teachings of the present invention. 

It is well-known that in x-ray CT fan-beam 
reconstruction certain factors can be applied to projection 
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data to diminish the effects of patient motion. While the 
mathmatics suggests that the present invention will not work 
with such prior methods, experimental results have 
demonstrated that some improvement is in fact obtained when 
the present invention is used in combination with such 
technique . 
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1. A method for producing an image, of a patient, the steps 
comprising: 

performing a scan in which a plurality of sets of 
projection data are acquired, each set of projection 
data being a view of a physical characteristic of the 
patient as seen from a different projection angle (9) ; 

acquiring a motion parameter with each set of 
acquired projection data, each acquired motion parameter 
being indicative of the position of the patient as the 
projection data is acquired; and 

back projecting each set of projection data to 
produce the image of the patient using a back projection 
formula which is modified by a displacement factor that 
has a value which is a function of the measured motion 
parameter; 

whereby artifacts produced in the image by movement 
of the patient from view-to-view during the scan are 
reduced, 

2. The method as recited in claim 1 in which the motion 
parameter acquired with each set of projection data 
indicates the position of the anterior chest wall of the 
patient. 
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The method as recited in claim 1 which includes: 

calculating a weighting factor for each set of 
projection data which is a function of the motion 
parameter acquired with the projection data and the 
projection angle (9) ; and 

correcting each set of projection data with its 
associated weighting factor prior to the back projection 
of the projection data. 

The method as recited in claim 3 in which the motion 
parameter acquired with each set of projection data 
indicates the position of the anterior chest wall of the 
patient and the weighting factor is calculated using the 
acquired motion parameter and a reference value which 
indicates the position of the anterior chest wall of the 
patient at a point in the patient's respiratory cycle. 

The method as recited in claim 1 in which the projection 
data is x-ray attenuation values. 

The method as recited in claim 5 in which the view-to- 
view movement of the patient is due to respiration and 
the weighting factor is calculated based on a geometric 
model of the motion of the patient's chest cavity during 
respiration and the acquired motion parameter which 
indicates the position of the patient's anterior chest 
wall. 

The method as recited in claim 6 in which the geometric 
model of patient motion is given by the formula 

f f U,y) - f(CX x + P x x, Oy + Pyy) 

where ji* and Py are magnification factors and a* and Oy 
are shift factors which indicate the point about which 
magnification occurs. 
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8. The method as recited claim 1 in which the set of 
projection data is acquired with parallel x-ray beams. 

9. The method as recited in claim 1 in which the set of 
projection data is acquired with a curved array of x-ray 
detectors that receive a fan-beam of x-rays, 

10. The method as recited in claim 1 in which the set of 
projection data is acquired with a flat array of x-ray 
detectors that receive a fan-beam of x-rays. 

11. Apparatus for producing an image of a patient, 
comprising: 

means for performing a scan in which a plurality 
of sets of projection data are acquired, each set of 
projection data being a view of a physical 
characteristic of the patient as seen from a 
different projection angle (9) ; 

means for acquiring a motion parameter with each 
set of acquired projection data, each acquired motion 
parameter being indicative of the position of the 
patient as the projection data is acquired; and 

means for back projecting each set of projection 
data to produce the image of the patient using a back 
projection formula which is modified by a 
displacement factor that has a value which is a 
function of the measured motion parameter; 

whereby artifacts produced in the image by 
movement of the patient from a view-to-view during 
the scan are reduced. 

12. Apparatus for producing an image substantially as hereinbefore 
described with reference to the accompanying drawings. 
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A method for producing an image substantially as hereinbefore 
described with reference to the accompanying drawings. 



